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Abstract 

Within the context of non-Hermitian quantum mechanics, we use the generators of eigenvectors 
of the Hamiltonian to construct a unitary inner product space. Such models have been of interest 
in recent years, for instance, in the context of VT symmetry, although our construction extends 
to the larger class of so-called pseudo-Hermitian Operators. We provide a detailed example to 
illustrate the concept and compare with known results. 
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Let us consider a non-Hermitian quantum mechanical Hamiltonian operator, H, which is 
related to its adjoint through a similarity transformation 

H = S'^H^S. (1) 

An operator with the above relation is sometimes called pseudo-Hermitian The adjoint, 
H\ is defined with respect to an auxilary Hilbert space, Ti, equipped with an inner product 
(■|-) (conventionally taken to be the Dirac inner product), such that {(plHip) = {H^^ (plip) , 
with lip), 10) G 7i. Relation ([T]) implies that 5* can be chosen to be self-adjoint. 

Since H is not Hermitian, it follows that e~**^ (for t G M, = 1) is not unitary and, 
therefore, it is not possible to construct a unitary quantum theory on 7i. To construct a 
unitary quantum theory we must define a new Hilbert space, 7i', with a modified inner 
product leading to a modified adjoint, H^, for which the Hamiltonian H is self-adjoint. 
By choosing Ti' so that = H, we can guarantee that e~**^ is a unitary time evolution 
operator in H'. 

Before we begin, let us make the following observation on the quadratic form in Ti defined 

as 

{<P\^)s := {<l>\S\ij). (2) 
The adjoint with respect to this quadratic form is given by 

{cP\H^l;)s = {H*ms. (3) 



where, 



It follows from ([H) that 



H* = S~^H^S. (4) 



H* = H, (5) 



so that H is self-adjoint, and time evolution is unitary with respect to this quadratic form. 
However, if the eigenstates of the Hamiltonian are given by 

= E\^,), Ee€, (6) 

then, it follows that 

{E - E'){tP,,\^,)s = {^,'\S{H - H)\'4j,) = 0, (7) 



where E denotes the complex conjugate of E. Equation ([7]) imphes that 

(V',^,), = e-^^(^)5,^., (8) 

where the constant e~^^^^'^ depends only on E because of the delta function, and may be 
negative (for instance, if 5* has negative eigenvalues). Consequently, while the quadratic 
form ([2]) leads to a unitary time evolution, it cannot be considered an inner product. 

We can turn this quadratic form into an inner product if we can remove the phase e~^^^^^ 
in (IHl), which will ensure that (V'b|V'£;)s > for all \iIje) G 'H'. If S* is a positive operator, 
so that it can be written as 5* = g^'g, then the quadratic form ([2]) is positive definite and, 
therefore, an inner product. (This will be the case, for example, when H is related to a 
Hermitian Hamiltonian h through a similarity transformation H = g~^hg, h'^ = h.) We will 
consider the case where 5* is not positive. We propose that there exists an operator A, which 
commutes with H such that 

q = SA, (9) 
is positive on a suitably defined Hilbert space, H', 

m), ■■= = {4>\SA\^) > 0, for all |0) G H'. (10) 

Since [A, H] = 0, we have 

:= q-^H^q = A'^S-^H^SA = A-^HA = H, (11) 

where denotes the adjoint with resect to {■\-)q- 
We define the spectrum of H to be 

spect(//) = {E s.t. < \ {^e\^e)s\ < oo}, (12) 

and note that ([1]) implies H and W are isosepctral. Let the projection operator P^, for 
E G spect(if), satisfy the following identities 

PeI^e') = 0, E^E', (13) 

and furthermore, let 

E£ spect(_H') 
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We will define Ti' =Im P^/ (the image of Pn') with, 

In' = Pn'- (15) 

Ti' equipped with {-\-)q becomes an inner product space once we define A which makes q 
positive. It follows from ffTTj) that e~**^ is also unitary for real t. 

We note that by construction {I^/'b)} is complete in H' and, therefore, any operator A 
that commutes with H can be expressed as 

A = Y,CePe, (16) 

E 

where G C, and it is understood that E G spect(if). It follows that 

q = Sj2cEPE, (17) 

E 

which, upon using (IHl), leads to, 

(^b'IV'e)? = {iJE'lS^CE" = CeO^e'ISHje) = Ce e''^^^'^ See' - (18) 

e" 

In order to ensure that the right hand side is positive for all {iPe) G 'H' we choose 



Ce = e 



iF{E) 



{'^E\ijE)-s\ (19) 



thereby defining the action of A. It follows that 

E B 

which gives 

{'ipE'li^E), = See', for all 1^^) G n'. (21) 

Comparing ([9]) with (120|) (and using (ITS!) ). we determine that A = e*^^^-* = (V'B|V'i5)s^ 
[E ^ H with if on the right), whenever {iPe\'4^e)s is a smooth function of E. However, the 
substitution E H may not be straightforward to carry out in practice, perhaps due to 
discontinuities in {ipE\'ipE)s (or as we shall see in the example, if the Hamiltonian describes 
more than one particle). For this reason we describe an alternative method of solving for q, 
assuming that one has already determined the operator which generates the eigenvector 
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Let us consider the following operator equation which defines as- 

Ha^ = Ea^ + a^kE. (22) 

We shall call a "generator" for the eigenstates of H if the following three conditions are 
satisfied^*! . 

{{) There exists at least one vector lip), solving 

k^\iP) = 0, for all E G spect(if), (23) 

with asl'ip) 7^ 0. 
(ii) There exists at least one vector |0) solving 

4 10) = 0, for all E E spect(i7), (24) 
with al~^\(j)) ^ 0, and {ip\(f)) ^ 0. 

{Hi) has an inverse a~^, at least on the subspace Pg, and an adjoint-inverse cr]j~^, at 
least on 16) defined in fl24l). 



For instance, if one has solved for r which satisfies the commutation relation [r, H] = i, 
then = e*^'^ satisfies the three conditions with = {H — Al), corresponding to the 
reference state lip) = {ipx), the state with energy eigenvalue A. In general, we can expect cTb 
to converge only for values oi E E spect(if) |3|. 
Condition {i) implies that 

\^e)=oM. (25) 

is an eigenvector of H with eigenvalue E. Furthermore, taking the adjoint of ( |22l) and 
formally acting on both sides with cr];^^ we get 

ift 4-1 = ^4-1 + 4-1 4. (26) 
It follows from condition {ii) that 

I0.) = 4"'I0), (27) 



[*] If IV'e) is «-fold degenerate we require a a% for each degenerate state. 



is an eigenvector of with eigenvalue E. 

Finally, from fl25|) and fl27j) we note that the eigenvectors of with eigenvalue E can be 
written either as S\iI)e) = SaE\'^), or as Therefore, these expressions must differ by 

a multiplicative constant depending only on E, 

c',^|z^,) = 4-^10). (28) 

We propose that 

c'^ = ce'nE)^ (29) 

where c is a constant independent of E and which can be seen as follows. Indeed, acting on 
both sides of (1251) with (V'eI, and using ^ as well as we obtain 

4e-^(^) = (^^|4"V) = (^l44"V) = m)- (30) 
By conditions (i) and (ii), the right hand side is independent of E, so we get c'^ = 

We therefore define the action of q on the eigenstates of H as 

q\il)^) = ga^lV^) := 4^^I0), (31) 

This is our main result. If we let go : 1'^) 10) (for instance, go might equal S) and 
normalize ("010) = 1, then we can express g as, 

g = $^4"'goa^'P^, (32) 

E 

and this leads to 

It is worth emphasizing here that this construction of g applies equally well in the case when 
the energy eigenvalues are complex as when they are real. As in the previous case, we can 
make the replacement E ^ H (with the i^'s to the right) if it is possible to express g in a 
power series in E. This procedure is demonstrated in the following example. 

We note that the strongest restriction on cr^ comes from condition {ii). For instance, 
ii H = + V{x)^ where p is momentum, one might try cTe = '0'e(x)/-0(x), which gives 
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ks = f{E,x)p:^i^ where f{E,x) is a function depending on E and x. {ip{x) 



denotes a wavefunction.) Clearly, kEip{x) = for all E, but kl;(f){x) = implies that 
<f){x) = {{f{E,x))^^)*. We see that condition (ii) cannot be satisfied if we define cTb this 
way. (We can weaken condition (ii) so that we require instead {ip\(f){E)) > for all (j){E) 
with the property kl;(f){E) = and \ip) as in (z)-we will examine this as well as other issues in 
a future pulbication.) Furthermore, this method is quite handy if we know the exact energy 
eigenstates of the theory. However, there are only a handful of such soluble examples in 
physics. Most systems can only be solved perturbatively. In this case, the generator cr^ can 
only be constructed perturbatively to a given order from a knowledge of the perturbative 
energy eigenstates to that order. The example that we discuss below, namely, the Lee model 
^, can be solved exactly. However, in a later publication we intend to explore the solubility 
of in greater detail. 



Example: 







The Lee model corresponds to an interacting field theory describing the decay of a 
massive fermion to another (massive) fermion and a (massive) scalar where all particles 
are charge neutral. The running coupling constant, in this model, becomes imaginary at 
high energies 51 and as a result, this model has been studied recently 6| in the context of 
PT-symmetry 7|, |8|] with an imaginary coupling where the interaction is non-Hermitian. A 
quantum mechanical model is described by the Hamiltonian 

H = Ho + Hi = {mgO^e N^N)+ig{e^N^V + V^N9), (34) 

where V, N correspond to the annihilation operators for the two fermions while 6 represents 
the annihilation operator for the scalar particle (actually, all particles in this quantum 
mechanical model are odd under parity, see for the behavior of various operators under 
parity, V, and time reversal, T). The Hamiltonian flM|) is invariant under VT and satisfies 
the relation H = VH'^V, where V = V'' is the parity operator (P^ = 1) and in this case, we 
can identify S = V which is not positive. 

It is easy to see that the combinations of the number operators Ny + Nj^ as well as N^+Ng 
commute with the Hamiltonian fl5^ . As a result, the energy eigenstates can be classified by 
the quantum numbers Uy + = 0, 1, 2 (as well as by n + = n,n + 1 where n denotes 
the eigenvalue of Ne). The states with Uy + = 0,2 are annihilated by the interaction 
Hamiltonian Hj (they correspond to trvial eigenstates of Hq) and, therefore, the nontrivial 



dynamics is contained in the eigenstates with ny + n^v = 1. For any finite n, this defines 
a two dimensional subspace spanned by the following two (unnormalized) eigenstates of H 
(states are labelled as |n, nv^, njy)): 

I^^J = a\n, 1, 0) + P\n + 1, 0, 1) = [a {O^YV^ + /3 (^t)n+iivt) |o), 
\^eO = {n + l)P |n,l,0) + a |n + 1,0,1) 

= {{n + l)P (^t)nyt + ^ (^t)n+i]vt) |0), (35) 

with energy eigenvalues 

= ^ ((2n + l)m0 + mN + mv- , 
K = \ ((2^ + ^)mo + mN + mv+ ^fJF^A^^ . (36) 

Here a = ^"'"V^ ^'^'^^^ (3, (3 = i ~ (''^s + "^^v — i^v) and 

|0) denotes the vacuum state of the theory. We note that for small coupling, the energy 
eigenvalues are real. We also see from fl35|) that the generators for {"^En) and |^E'_b;) are given 
by 

aE„ = a {e^rV^ + P (^t)n+iArt, 

aj,,^ = {n + 1)^ (^t)nyt + ^ (e^Y+^Nl (37) 

That is, = (Ti\0), for z = En, E'^. By evaluating [if, cr,;] we find ki oc H + Ng + Nj^, and 
therefore, fc||0) = fci|0) = 0. Thus, we can choose = \(f)) = |0) so that go = 1- 

Let us next solve for g„ which represents the action of q restricted to this two dimensional 
subspace. In constructing the inverse for aE„ we note that there is some freedom following 
from condition {iii), which only requires the inverse to be well defined on P„. We use this 
freedom to make the additional requirement that a^^l'^E^) =0. In this way we avoid 
constructing a projection onto l^^En)- Doing the same for cr^/^ we obtain, 

(TeI = -4 + ^^TTT (38) 
n\ (n + l)\ ^ ' 
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It can be checked that — (^ij|0), as claimed. Let us look at the product 

n\ J 



ni \ ni 



TV. \n\ n\ J 



= — [\a\^ Ny{l - N^) + (n + l)iV^(l - Ny) - a(3 ONV^ - a(3 O'^N'^V] . (39) 



The last relation follows from the fact that the expression acts only on the states |n, 1,0) 
and \n + 1, 0, 1), where, for instance, we can use {9^Y6'^\n, •, •) = n\\n, •, •). The constraints 
Nv + Nr, = l (in this space) as well as A^^ = = Nl, allow us to replace F V by Ny{l-N^) 
and N'^N by iV^(l - A^v-)- Likewise, obtain 

-a^ 9NV^ - a(5 O^NW] . (40) 

We can sum these two terms, and then make the replacement n — >■ O'^O = Ng to obtain 
q. Before doing so, we will exploit the freedom to multiply each g„ by an arbitrary positive 
constant. Doing this will change the normalization (^j|^j)^ = 1, but the innerproduct will 
remain positive. Therefore, instead of defining g„ = {^eIV^eI + (^e^V^e^^ simplify 
the expression by multiplying out the factorials dividing the previous two equations (giving 
instead {^f eJ^/ eJ ^ = n\) . 

With this in mind, we use the identities jap + (n + l)|/5p = , ^ ^=, and /3a = 



, ^=, to obtain 

VM2-4fl2(n+l)' 



nl i^Ey^El + in + l)H^Eiy^El 

(|q;|2 + {n + 1)|/3|2) {Ny{l - N^) + N^{1 - Ny)) - 2ap ONV^ - 2(3a O^N^V 

^ {Ny{l-N^)+N4l-Ny)) 



^/fi^-4g\n+l) 
+ , _ ^^^^ _ 2zg 

In order to make the expression valid for an arbitrary number of scalars we replace n by 
^t^ — We must be careful about the operator ordering, for instance, 9^ increases 

the scalar particle number by one. In addition, the fermion operators V, N annihilate the 



groundstate so that we should add a projection onto |0) given by Pq = {1 — Npf){l — Ny) 
to finally get: 



q = Po + J2^^ — >in^ Ne) 



+eNV^ , - , e^N^V. (42) 

V/i^ - 4(7^ AT, _ 4^2jv^ 



The above expression was previously determined [6| using a more elaborate method. In 
that case, the result was derived, in part, by requiring that the operator C = qV is an 
involution (C^ = 1) when acting on the eigenstates of H (it is actually a projection onto the 
set {|vl>sj, |vl>s;)}). 

In summary, we have utilized the fact that the generators of eigenstates of a non- 
Hermitian operator, H, contain sufficient information for constructing a positive inner prod- 
uct space Ti' . We have illustrated how the method works in the case of the VT symmetric 



Lee model 



6|. In a subsequent publication we plan to consider possible constructions for 



the operator cr^;, as well as its role in more sophisticated models. 
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